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Abstract 

We study the full unitary matrix models. Introducing a new term llogU , I 
plays the role of the discrete time. On the other hand, the full unitary matrix 
model contains a topological term. In the continuous limit it gives rise to a 
phase transition at 6 = tt. The ground state is characterize by the discrete time 
I. The discrete time I plays like the instanton number. 



1 



1 Introduction 



Models of the symmetric unitary matrix model are solved exactly in the double 
scaling limit, using orthogonal polynomials on a circle. [0 The partition function 
is the form J G?f/exp{ — j-tiV(U)}, where U is an N x N unitary matrix and 
trV(U) is some well defined function of U. When V(U) is the self adjoint we 
call the model symmetrical The simplest case is given by V(U) = U + . This 
unitary models has been studied in connection with the large- N approximation 
to QCD in two dimensions. (one-plaquette modet)^ The full non-reduced uni- 
tary model was first discussed in ||. The simplest potential is divided into the 
symmetric and anti-symmetric part.Q The symmetric part is the usual Wil- 
son action. The anti-symmetric part becomes the topological term.|^| In the 
continuous limit we can reconstruct the topological information encoded in the 
theta term. We introduce a new term I log U. The full unitary matrix model 
can be embedded in the two-dimensional Toda lattice hierarchy. Q Adding the 
new term to the full unitary matrix model, I becomes the discrete time and 
the partition function satisfies the discrete time Toda equation. This equation 
is sometimes called as "Hirota equation" . JtJ Discrete time integrable system is 
currently under extensive study. The partition function for the Hermite matrix 
model with the new term also satisfy the discrete time Toda equation. || In 
the full unitary matrix model the discrete time I couples to the theta term and 
characterizes the ground state. The partition function for the SU(N) case is 
the sum of the partition function for the U(N) case about all the discrete time 
I. Thus the theta term would have been zero for SU(N) case. 

This letter is organized as follows. In the section 2 we introduce a new term 
to the full unitary matrix model. We show that the partition function of our 
model satisfy the discrete time Toda equation. In the section 3 we see the role 
of the discrete time / in the Painleve equation. Coupling the Toda equation and 
the string equation, we obtain the Painleve V with 8 = 0. This becomes the 
Painleve III. In the section 4 we consider the physical meaning of our model. 
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We derive the continuum limit of the action which will turn out to be the heat 
kernel with the theta term. We interpret the continuum partition function in 
terms of the underlying topology and show the presence of a phase transition 
of at 8 — 7r. The ground state is indicated by the discrete time I. The latter 
half of this section we calculate the partition function for SU(N) and study an 
effect of the theta term. The last section is devoted to concluding remarks. 

2 The partition function on the lattice 

We consider the partition function of the full unitary matrix model with a 
discrete time. We calculate explicitly the determinant form of the partition 
function. 

We consider the unitary matrix model 

/N 

dUaq>(-- 5 V(U)), (2.1) 

where V(U) is a potential 

V(U) =t 1 U + t- 1 U* -llogU. (2.2) 
U is the gauge group U(N). Here we consider the case 

i-i 

We parameterize t\ and t-\ by e: 

ti = -e e , t-i = -e~ e . (2.4) 
The measure dU may be written as 

dU = f[^A(a)A(a). (2.5) 

m 

Here the eigenvalues of U are {exp(iai), exp(ia2), • • • , exp(iajv)} and AA is 
the Jacobian for the change of variables, 

A(a) = Act ]k e la ^ N - k \ 

A(a) = detj-fce - ^^-*). (2.6) 
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In terms of these variables the potential V(U) is 

N 

V(a) = — (cosh e cos ay + isinhesinctj +ilaj). (2-7) 

3=1 

Putting this all together and introducing 

f(a) — cosh e cos a + isinhesina + ila, (2.8) 



we obtain 



Z^N/g 2 ) = const. J] y^e^)|A(a)| 2 , 

N r* a N 
= const, n / ^e^lEMrn^ 



AT 

N 

const. 2 (-1)" E II 

AiGSiv cr£:Sn m=l ' 

x exp[/(a m ) - ia m (a(m) - /jcr(m))], 



2?r 



= const. AndewT — e /(««)-fa(*-i)]. (2.9) 

where the summations are on all the permutations of the N elements. We can 
calculate the (j, k) element of the determinant as 

r da < n , u . . i ■ \ • i 

/ — exp — (cosn e cos a + i smn e sin a) — ima\ 

= / — exp[^cos(a-ie)-ima] = e me I m (N/g 2 ), (2.10) 
where m = I — j + k. Here / m is the modified Bessel function of order m. 



Substituting this into (2.S) we can obtain the partition function 

Z^ x = const.AMdet^e^-^Vj+fcW^ 2 ). (2.11) 

Here we introduce symbols D, and D?'2' which are the determinants of 
the (N + 1) x (N + 1), N x N and (N - 1) x (N - 1) matrices respectively with 
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the respective definitions. We define D = Z^ +1 ; . D£ is same as D except that 
the a-th row and 6-th column are removed from it. is same as D except 

that the a-th and c-th rows, and 6-th and d-th columns are removed from it. 
Then the Jacobi formula 

D\D N N X\ - D? +1 ■ Dh +i = D ■ D\:Zl (2.12) 
becomes the recurrent relations 

(^n,i) 2 ~ ^N,l+l ' %n,i-i — %N+i,l ' ^n-i,i (2-13) 

This is the discrete time Toda equation. I is the discrete time variable. In the 
continuous limit of I, (2.13) becomes the Toda molecule equation. ( 2.13| ) is the 



same relations which is obtained from the symmetric model. |j) 

3 Unitary Matrix Model and Painleve III 

It is well known that the partition function Zjlj of the unitary matrix model 
can be presented as a product of norms of the biorthogonal polynomial system. 
Namely, let us introduce a scalar product of the form 

<A : B>=j ^1 exp{-V(z)}A(z)B(z- 1 ), (3.1) 

where 

V(z) = (tmz m + t. m z- m ) - I log z. (3.2) 

m>0 

Let us define the system of the polynomials biorthogaonal with respect to this 
scalar product 

<<f>n,h$t,t>=h k ,i5 nk . (3.3) 
Then, the partition function Zj^ ( is equal to the product of /i n 's: 

JV-l 

z u n.i = n h ^> t ° =i - ( 3 - 4 ) 

fc=0 
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The polynomials are normalized as follows (we should stress that superscript '*' 
does not mean the complex conjugation): 

$n,i = z n + ■ ■ ■ + S n - hl , K,l = z n + ■ ■ ■ + S- U = S*_ ltl = 1. (3.5) 

Now it is easy to show that these polynomials satisfy the following recurrent 
relations, 

$ n+ l,;(z) = Z$ n .l{z) +S n .lZ n $* n j{z- 1 ), 

K+iAz' 1 ) = ^"^^(^l + ^-AlW, (3.6) 

and 

(3-7) 

Note that h n j, S n j, 5 1 * z ,<& n j(z) and , depend parametrically on t±,t2, 
and i_i,i_2, • • • , but for convenience of notation we suppress this dependence. 
Hereafter we call ti, t%, • • • , and t—i, i-2, • • • 5 time variables. 

Using ([si]) and integration by parts, we can obtain next relations: 



27riz 

and 

d/i(z 



V'(z)9 n+1 , l {z)9* l (z- 1 ) = (n+ l)(h n+u - h n ,i), (3.8) 



2mz 



z 2 V'{z)$* n+1 iiz- 1 )®^) = (n + l){h n+u - ft,,,,). (3.9) 



(p.q) and (3.9) are string equations of the full unitary matrix model. 



If t\ and t-i are free variables while ti, = is = ■ ■ ■ = and t-2 = £-3 = 



0, (3.5) and (3.9) become 



(3.10) 

(n+l)S n ,lS* t i = ti(S* l S n +ld+SndS* l _ ll )(l — Sn i lS* l )+l(l — Sn,lS^ ll ). (3.11) 



Next we introduce a useful relation. Using (3.3) and integration by parts, 
we can show 

d(i(z) 



2-kiz 



■zV'(z)<P n j(z)<S>* l (z- 1 ) = 0. (3.12) 



This corresponds to the Virasoro constraint: 







^ ot r , 



(3.13) 



k— — oo 

This relation constrains a symmetry like the complex conjugate between and 
i_fc. If we set that t\ and i_i are free variables while t% = i3 = • •• = and 



(3.14) 



t_2 = t-3 = • • • = 0, from ( J3 . 1 2| ) we get 

tlSn,lSn~i,l = t-iS^iSn-iJ — I. 

Using (|TlJ), ( |3.10| ) and fl3.1l| ) can be written 

(n + l)S„ t i = (t-xS n +i,i + txS n -x,i)(l - S n< iS* t i), (3.15) 

(n + 1)5:, = + *-iS:-u)(l - Wn,«)- (3.16) 

Using the orthogonal conditions, it is also possible to obtain the equations 

which describe the time dependence of Q n ,i(z) and Namely, differ- 



entiating (3.3) with respect to times t\ and i_i gives the following evolution 
equations: 



dtx S n -x,l h n -x,i 



dt-x h n -x,i 
SKA*- 1 ) h. 



$n-u( 2 )> 



dtx 



s 



1,1 



-($* Az- 1 ) - z-^.xj), 



(3.17) 
(3.18) 
(3.19) 
(3.20) 



dt-x S n-x,l K-l,l 

The compatibility condition gives the following nonlinear evolution equations 



dS n ,i c h n +i,i 
— —&n+i,r 



dtx 
dS* 



h n ,l 

h n +i,i 



9S n> i h n+1 ,i 
— -iti-i,r 



dt-x 

dS h__ q , 



h n ,i 

h n +x,i 
h n ,i 



~~dtT ~ s:+u 

9h n j dh n j 

— <-'ra,Z>-' n _l i>l"n,,li "TT, — ^ n l&n — 1,1 "n ,1 ? 



dtx 



dt-x 



(3.21) 

(3.22) 
(3.23) 
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Here we define <x n j, b n ,i and b 



a n,l = 1 — S n> lS n i 



ln+l,l 



(3.24) 

(3.25) 
(3.26) 



b n ,i — S n jS n _ ll , 

= S nt lSn-l,l- 

Notice that since the definitions o, n b n i and 6* ; satisfy the following identity: 

»M n,i = (1 - - On-1,0- (3.27) 



It can be shown using ( 3.14|) that 



(3.28) 



In terms of a n> i, b n< i and b* n t , ( 3.21 ) and ( 3. 22] ) become the two-dimensional 
Toda equations: 



da n j db n ,i 

a>n,i\O n +i t i — Onj), — — = a n ^i — a n -i,i, 



and 



dti 

da„j 
dt-i 



a n ,i(b 



n+l,l u n,U> 



dt-1 

dtl 



Using a n j, b n ,i and b* nl , we rewrite ( 3.10 ) and ( 3.11 ) 
n + 1 1 - a n j , , , 



and 



n + 1 1 - a n j 



t-\ 0>n,l 

From (U§ and dlj| ) we eliminate o n +i .;, 

26 n ,z = [— (1 - On,0 - ^— ] 



In the same way, from ( 3.30 ) and ( 3.32 ) we eliminate b 



n+l,/' 



2b n,l = [-7 (1 - On,;) " "77T J- 



(3.29) 



(3.30) 



(3.31) 



(3.32) 



(3.33) 



(3.34) 
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Using ( p^27| ) and ( ^28| ), fT^ ) and pT30| ) can be written 



db, h i 1 b n .i(tib, h i +1) 

- (a n ,i - 1) + 



56: 



9ti 



= (a„,i - 1) + 



t-i 1 — a n ,i 

1 b^Ct-lftn.J - 



*1 



1 - ou 



(3.35) 



(3.36) 



Using (3.33) and (3.35) to eliminate we obtain a second order ODE for a„ 



d 2 a ri 



n + 1 <9et„ i n + 1 <9a,u 



2a„,i(a nj i - 1) + 



£-i»n,2 9ti £ia n ,z 9t-i 
a n ,i dti dt-i 2 t_i (a n j - l)a n ,i dh 



(n+ l) 2 a n j - 1 

wi — ^f) 

Lt\t-\ a n .i — 1 



In the same way, we eliminate 6* ; using ( 3.34 ) and ( 3.35 ) and obtain an ODE 
for a„i 



d 2 a n 



dhdt^ 



n+1 da U: i _ n+l da n ,i 
t\(in,i dt-% t^ia n j dti 
1 da n i da n , , 1 t_i 



2a, h i(a n ,i - 1) + 
1 



(n+ l) 2 q„j - 1 

; 2 



Using (|H|), ( |3.29| ) and fl3.30| ) we can obtain 



^3§) 



h ~mT = 

So a„ / are functions of the radial coordinate 



(3.39) 



x = tit 



lE-lj 



(3.40) 



only. Then from ( 3.37 ) and ( 3.38| ) we can obtain 



d 2 a n ,i _ 1 , 1 

dx 2 * [ 



2 a n ,i - 1 a n ,i dx 



1 w 9a«j s 2 1 d fl ru 



x dx 



(n + l) 2 a,, ht - 1 f 



--a n ,l{ a n,l - 1) + 9 2 



2x 2 a r , 



1 



(3.41) 



( 3.41 ) is an expression of the fifth Painleve equation (PV) with Sv = It is 
interesting that (3.41) has a symmetry a n j <-* a„.z — 1 and (n + l) 2 <-> — Z 2 . To 
see this, we define c n X- 

c n ,l = l-a«J- (3-42) 
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Using Cnj, through the transformation x — ► — x we rewrite ( 3.4l| ) 

d 2 C n j _ 1 1 1 w gCn£ v2 1 dc ra ,; 

9x 2 2 c n ,i — 1 c Ut i dx x dx 

1 Cnl ( Cnl n i - c ">' -1 (" + 1 ) 2 c ",' ( 3 43 ) 
x 2x 2 c n i 2x 2 / 1 



In (3.43) the role of I 2 and (n + l) 2 exchange. 

This equation can be transformed into the third Painleve equation (P III), 
(see 0,0]) ( 3.41 ) is the same equation obtained from the 0(4) nonlinear a- 



model.[|l). 

4 Topology on the lattice 

We divide the potential into the symmetric and the anti-symmetric part. 

V(U)=t+s w +t^s e , (4.1) 

where 

s w is the symmetric part, the usual Wilson action 

s w = htrU + trtf). (4.3) 



Here we choose 



8fl = i(trl7-trl7t), (4.4) 



for the theta term, the anti-symmetric part. In the naive continuum limit, the 
sum of the theta term gives the first Chern number of the bundle on which the 
gauge field lives. We parametrise the one plaquette action as 

s(U) = -coshes w {U) - smhes e {U) - llogU, (4.5) 

where e is a real parameter that determines the relative weight of the two terms, 
introducing in the section 2. 
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We calculated the partition function in (2.11). We consider the continuum 
limit of the partition function to leading order in g — > 0. As g — » we want the 
coupling of the theta term to approach its continuum value 9/2ir. Accordingly 
e has to be scaled in this limit as 



2tt N' 

The asymptotic representation of the modified Bessel functions are 

6 Z n 2 

In{z) — > p~~ as z 



'2-kz 



oo. 



Substituting (4.6) for e and (4.7) into (2.11), we obtain the small g 



(4.6) 



(4.7) 



Zf{9,6) 



det jk exp[-f-{(Z -j + k) 2 --(l-j + k)} 
det jk e XP [-^{(/ - j + k f f -/-?}]. 



(4.8) 



Notice that 6* couples to the discrete time I. For our purpose we are interested in 
the coupling of 9 and I. Then in the first step we drop the independent factors. 
The action is given by 

8(0, l) = - log(det jfc exp[-^-{G -j + k-ff- -/-) 2 }]). (4.9) 

ZiV Z7T 4 7T 

If 9 is varied, the ground state might change at some points due to the level 
crossing which causes phase transition. 
For 9 = using next relation 



detjk exp —a(m — j + fc) 2 = e am N detjk exp —a(—j + fc) 2 , 



(4.10) 



the action is 

2 2 

s(0,l) = -log(det jk exp[-^(l-j+k)]) > -]og(det Jfc exp[-^(-j+A)]) = s(0. 

(4.11) 

It is clear that at 9 = the ground state is I = 0. 
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For any 9 using (4.10) we can obtain that 



•M«./) = -^(-) 2 -logdet jfc exp[-|^-{(/-j + fe--^ y 

7T ZiV Z7T 

J2 



8 + V ~ 2^ )2 " l0gdetjfe 6Xph lv ( ^ + fc )^ 4 ' 12 ) 

The first term is a constant shift of the free energy. The second term shows 
that by increasing 9 staring from zero we do not encounter any phase transition 
up to 9 = 7r. At this point however there is a phase transition. If 9 is above 
7r, I = 1 becomes the ground state. Due to the periodicity in 9 the same type 
of phase transition occurs whenever 9 becomes an odd integer times 7r.(Fig.l.) 



From (4.12) the free energy density around the phase transition is given by 



f(0) = 0, -7T>9>7T, 
2 

f(9) = *L(l-~) - 3tt > 9 > -7T, 3tt > 9 > tt, 
2 p 



7T 



/(0) = g 2 {2--) - 5tt > 6» > -3tt, 5tt > (9 > 3tt, (4.13) 



We have to note that ( 4.13 ) gives only the part of the energy without the 
constant term. 

Next we consider the SU{N) case. In this case there is an additional 
X}fc=i a k — constraint on the eigenvalues that must be enforced by a delta 
function in the measure. If the delta function is written as 

N oo N 

fc=l m=-oo fe=l 



The calculation (2.E) goes through without any modification giving the final 
result, 

oo 

Zffi = const.AH J] det jk e^ l+m -^I l+m ^ +k (N/g 2 ) 

771 — — OO 

OO 

= const.AH det jk e e< - m ^ +k h m ^ +k {N/g 2 ). (4.15) 

m— — oc 

In the last step the discrete time / disappears for the summing about m. Using 



12 



(4.7) for small g 1 we obtain 

oo oo 

Zi U (9,0) = E const.AT! £ det jk e< m - j+ ^I l+m ^ +k (N/g 2 ) 

771 — — OO 771 — — OO 

~ E det jfe exph^{(m- J + fc--) 2 - i (-) 2 }] 

m= — oo 

00 2 n 2 ia 

= Y, ex P[-y( m -^) 2 ] det ^ ex P[-|^{(-J + fc ) 2 -i(^I 

m— — oo 

In the g — > limit this sum can be explicitly evaluated by the approximating it 
with a Gaussian integral. Then we can obtain 



Z? u (g,9) ~ ^_exp^(-) 2 det, fe exp[-^{( J -fc) 2 - -(-) 2 }]. (4.17) 
Then in the 5/7 (iV) case has influence on a constant shift of the free energy 



only. We can rewritten (4.15) as 

oo 

Zf% = const.TV! £ tet jk e< l - i+ie > b-j+^N/ ff 1 ) 

I — — OO 

oo 

= E ( 4 - 18 ) 

/ — — oo 

The partition function for the SU(N) case is the sum of the partition function 
for the U(N) case about all the discrete time I. Thus the theta term would have 
been zero for SU(N) case. 

5 Concluding remarks 

In this letter we introduce a new term to the full unitary matrix model. This 
term is I log U. We show that the partition function of this model satisfy the 
discrete time Toda equation and I is the discrete time. Furthermore we consider 
the relation the partition function of this model and Painleve III. We can see 
the symmetry between n and /. 

We derive the continuum limit of this model in the presence of a topological 
term. The theta term remained to be coupled to the dynamics in the continuum, 
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giving rise to phase transition at 8 C = ±7r, ±37r, ±57r, • • •. The ground state is 
indicated by the discrete time I. The discrete time plays like the instanton 
number. The partition function for the SU (N) case is the sum of the partition 
function for the U(N) case about all the discrete time I. Then for the SU(N) 
case the theta term has influence on a constant shift of the free energy only. 
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Figure 1: The free energy density around the phase transition. 
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